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Studies of human germ cell development are limited in large part by inaccessibility of germ cells during devel- 
opment. Moreover, although several studies have reported differentiation of mouse and human germ cells 
from pluripotent stem cells (PSCs) in vitro, differentiation of human germ cells from PSCs in vivo has not 
been reported. Here, we tested whether mRNA reprogramming in combination with xeno-transplantation may 
provide a viable system to probe the genetics of human germ cell development via use of induced pluripotent 
stem cells (iPSCs). For this purpose, we derived integration-free iPSCs via mRNA-based reprogramming with 
OCT3/4, S0X2, KLF4 and cMYC alone (OSKM) or in combination with the germ cell-specific mRNA, VASA 
(OSKMV). All iPSC lines met classic criteria of pluripotency. Moreover, global gene expression profiling did 
not distinguish large differences between undifferentiated OSKM and OSKMV iPSCs; however, some differ- 
ences were observed in expression of pluripotency factors and germ cell-specific genes, and in epigenetic pro- 
files and in vitro differentiation studies. In contrast, transplantation of undifferentiated iPSCs directly into the 
seminiferous tubules of germ cell-depleted immunodeficient mice revealed divergent fates of iPSCs produced 
with different factors. Transplantation resulted in morphologically and immunohistochemically recognizable 
germ cells in vivo, particularly in the case of OSKMV cells. Significantly, OSKMV cells also did not form 
tumors while OSKM cells that remained outside the seminiferous tubule proliferated extensively and formed 
tumors. Results indicate that mRNA reprogramming in combination with transplantation may contribute to 
tools for genetic analysis of human germ cell development. 



INTRODUCTION 

A significant challenge in elucidating genetic requirements for 
human germ cell formation, maintenance and differentiation is 
to recapitulate germ cell specification and differentiation both 
in vitro and in vivo. Studies in the mouse suggest that complete 
reconstitution of mammalian germline development from pluri- 
potent stem cells (PSCs) is possible (1-4). In spite of successes 
in the mouse, in vitro differentiation of human PSCs to germ 
cells that progress through meiosis and form functional gametes 
remains a significant challenge. Notably, previous efforts, 



including our own, have used a variety of in vitro methodologies 
that advanced studies of human germ cell differentiation but 
consistently yielded low numbers of germ cells, inconsistency 
across line derivations and genotypes and incomplete imprint 
erasure and re-establishment in a sex-specific manner. Here, we 
sought to differentiate human germ cells by directly transplanting 
undifferentiated human induced pluripotent stem cells (iPSCs) 
into murine seminiferous tubules in order to make use of the 
germ cell niche to promote human germline formation in vivo. 

Insights gained from mouse studies, combined with knowl- 
edge of genetic control of PGC fate in diverse organisms and 
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our previous studies of human germ cell differentiation in vitro, 
prompted us to consider whether developmental 'priming' by 
induced expression of germ cell-specific genes in combination 
with xenotransplantation might overcome current roadblocks 
to efficient human germ cell development from PSCs. In 
murine studies, a critical advance was made by the induction 
of a transient epiblast-like cell state from mouse embryonic 
stem cells (mESCs) and iPSCs by addition of cytokines or via 
overexpression of PRDM14, DPP A3 and TFAP2C transcription 
factors alone (2-4). Resulting cells were developmentally com- 
petent to form spenn and oocytes following transplantation to 
produce live, healthy offspring after fertilization. In our own 
efforts, we have used overexpression of key germ cell regulators 
including DAZ gene family members or DAZL and VASA to- 
gether to drive germ cell differentiation and meiotic progression 
from human ESCs and iPSCs in vitro (5-7). However, our 
studies and those of others using in vitro mediated differentiation 
have been confounded by low yields of germ cells, inefficient 
meiotic progression and an incomplete imprinting status (8,9). 
Owing to inherent differences between human and mouse PSC 
and, based on previous studies, we predicted that induced expres- 
sion of translational regulators such as VASA, DAZ, DAZL and 
BOULE might promote human germ cell formation. Thus, we 
included VASA, a translational regulator, to the mix of factors 
used in mRNA reprogramming to iPSCs in hopes of alleviating 
hurdles that we and others have encountered with human germ 
cell derivation in vitro (5-9). The VASA gene encodes a highly 
conserved germ cell-specific RNA-binding protein whose role 
in germ cell development may include acting as a chaperone to 
enable conxct folding of different target RNAs in germ cells 
(10). Furthermore, we note that similarities between pluripotent 
human ESCs and iPSCs to mouse epiblast cells lends support to 
our rationale that we might produce primed iPSCs for germ cell 
development (11-14). We then transplanted the undifferenti- 
ated iPSCs directly into the seminiferous tubules of germ cell- 
depleted immunodeficient mice, in order to evaluate the contri- 
bution of VAS A-primed and non-primed cells to germline devel- 
opment in vivo. 

RESULTS 

Addition of VASA to OSKM-reprogramming mix 
does not affect reprogramming efficiencies or kinetics 

Synthesized modified mRNAs, including VASA mRNA, were 
produced and validated prior to use as shown (Supplementary 
Material, Fig. SI A and B); protein expression and localization 
was confirmed by immunocytochemistry (Fig. 1). We added 
synthesized VASA mRNA to OSKM factors in molar ratios of 
3:0.5:1:0.5:1 (OSKMV) and examined reprogramming effi- 
ciency and kinetics of different somatic lines [BJ (XY), HUFl 
(XY), HUF3 (XX) and HUF9 (XX)] with either the OSKM or 
OSKMV cocktail in parallel. We adapted a previously reported 
protocol (15) in which feeder- free iPSC derivation was accom- 
plished with six mRNA factors (OSKM NANOG 
LIN28A) and modified OCT3/4. We successfully derived 
feeder- and xeno-free iPSCs with OSKM alone and with 
OSKMV (Supplementary Material, Fig. S2A-F). We observed 
colony formation with both OSKM and OSKMV reprogram- 
ming with no significant differences in colony number 



(~0.5% efficiency) or timing of colony appearance (6-12 
days post first-transfection). We then further analyzed the XY 
lines, HUFl and BJ, and compared OSKM- and OSKMV- 
derived lines across a series of functional and molecular assays 
beginning with a comparison of morphology, alkaline phosphat- 
ase ( AP) staining and embryoid body (EB ) formation. We did not 
observe notable differences, regardless of reprogramming cock- 
tail used (Supplementary Material, Fig. S3A). Colonies repro- 
grammed with OSKMV were characterized by a high nucleus/ 
cytoplasm ratio, prominent nucleoli, well defined borders and 
a distinguishing chromatin structure and nuclear architecture 
(speckles and heterochromatin domains), all features that are 
very similar to OSKM-reprogrammed colonies and human em- 
bryonic stem cells (16). 



Transient ectopic VASA expression alters gene expression 
signatures of derived iPSC lines 

As part of the assessment of pluripotency, we examined en- 
dogenous gene and protein expression of various markers asso- 
ciated with pluripotency including P0U5F1, NANOG, SALL4 
and DNMT3B in both OSKM and OSKMV colonies (Supple- 
mentary Material, Fig. S3C). Notably, expression of a subset 
of markers associated with pluripotency was lower {P < 0.05) 
in lines reprogrammed with OSKMV relative to their OSKM 
counterpart, with PRMT5, SALL4 and DPPA4 being the most 
significantly different (P < 0.001). We also confirmed a 
similar reduction in expression of a subset of genes in lines 
that were derived with OSKM or OSKMV via a lentiviral repro- 
gramming strategy to exclude reprogramming strategy related 
events (Supplementary Material, Fig. S3C). We then examined 
effects of transient ectopic expression of VASA during repro- 
gramming on expression of genes associated with early germ 
cell development. We observed that the majority of markers 
showed gene expression levels similar to the lines 
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Figure 1. Functional validation of mRNA expression encoding for VASA. 
Immunostaining of VASA in mRNA-transfected fibroblasts 24 h after transfec- 
tion. VASA protein localized correctly in the cytoplasm. Mock transfected, sec- 
ondary antibody stained only and non-transfected samples served as negative 
controls. Scale bar, 10 |jLm. 
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reprogrammed with OSKM alone and/or the parental fibroblast 
line, indicating no gene activation (exemplified by PRDMl). 
However, a subset (PRDMl 4, DPP A3 [STELLA] and VASA) 
was expressed at significantly higher levels {P < 0.001) in 
iPSC lines reprogrammed with OSKMV relative to OSKM- 
derived colonies, indicated for iPSC.HUFl cells (Fig. 2A). 
Results were partially mirrored by the lentiviral-derived HUFl 
iPSC line and mRNA-derived iPSC.BJ cells. We detected a 
smaller subset of germ cell markers at a higher expression level 
in OSKMV-derived lines of which only PRDMl 4 (for iPSC.BJ) 
was significantly upregulated compared with the OSKM counter- 
part (Supplementary Material, Fig. S3D). We note that gene ex- 
pression was measured at two different passages (passage 4 and 
14) to eliminate the possibility of expression from exogenous 
mRNA and to demonstrate stability of the distinct endogenous 
gene expression profile. We further confirmed endogenous 
VASA gene expression in OSKMV cells with immunocytochem- 
istry (Supplementary Material, Fig. S3B). 

OSKMV-derived lines, similar to OSKM-derived lines, 
are fully pluripotent 

To further assess whether mRNA-reprogrammed OSKM and 
OSKMV lines are both fully pluripotent, we performed a 
variety of molecular and functional assays. Spontaneous dif- 
ferentiation revealed that all lines form all three germ layers 
in vitro (Supplementary Material, Fig. S3E). Epigenetic analysis 
revealed similar methylation status of OCT3/4 and NANOG 
promoters of both OSKM- and OSKMV-derived HUFl lines 
comparable to hESCs, which were substantially different rela- 
tive to the parental fibroblast line (Supplementary Material, 
Fig. S3F). Moreover, we observed that all clones were karyoty- 
pically normal and differentiated to all three germ layers in vivo 
in teratoma assays, following kidney capsule injection into SCID 
mice (Supplementary Material, Fig. S3G and H). Taken to- 
gether, these results indicate that mRNA-derived iPSC clones 
from multiple independent fibroblast lines, with and without 
addition of VASA are fully reprogrammed by the most stringent 
assays available to date for human pluripotent cells (17). 

OSKMV- and OSKM-derived clones respond differently 
to BMP4 treatment 

VASA overexpression from lentiviral vectors has been reported 
to drive differentiation into primordial germ cells (PGCs) and 
more mature germ cells when overexpressed in both human 
ESCs and iPSCs (6). In mice, PGC specification is characterized 
by localized BLIMP 1 (PRDM 1 ) expression in a small number of 
cells in the posterior, proximal mouse epiblast at E6.25 followed 
by upregulation of DPP A3 around E7.25 (18). Several groups 
have demonstrated in vitro differentiation of PGCs from both 
mouse ESCs and hESCs via media supplementation with bone 
morphogenetic proteins (BMPs) (3,5,7,19). Thus, we hypothe- 
sized that lines derived with OSKMV might reveal a different 
potential in their ability to respond to BMP4 relative to OSKM 
cells. We supplemented culture media of both OSKM- and 
OSKMV-derived lines with BMP4 over the course of 4 days 
and assessed activation of gene expression of early germline 
markers. We included two controls — an mRNA-derived 
(OSKM) BJ fibroblast line and a hESC line (H9). Although a 



subset of genes, including NAN0S2, RET and Y chromosome 
DAZ2 among others, did not show significant changes in gene 
expression in either the OSKM- or OSKMV-derived lines, 
expression of NAN0S3, VASA and DPPA3 was significantly 
upregulated (P < 0.01) post-BMP4 treatment in the line derived 
with OSKMV after 4 days of BMP4 treatment (Fig. 2B). iPSCs 
derived with OSKM as well as controls did not demonstrate 
differential expression of these germ cell markers in response to 
BMP4 treatment, indicating that observed differences are con- 
fined to different reprogramming factor cocktails used during 
iPSC derivation. 

We also examined gene expression changes of pluripotency- 
associated markers during BMP4 supplementation (Supplemen- 
tary Material, Fig. S4A and B). DNMT3B, P0UF5 1 and SALL4, 
all genes that are strongly expressed in hESCs and iPSCs, showed 
a rapid decrease in gene expression after 4 days of BMP4 treat- 
ment in both OSKM- and OSBCMV-derived lines as well as both 
control lines, as expected (Supplementary Material, Fig. S4B). 

OSKMV- and OSKM-derived lines display distinct 
SCP3 staining patterns 

Synaptonemal complex protein 3 (SCP3) encodes a meiosis- 
specific protein that is essential for formation of meiotic syna- 
ptonemal complexes of the maternal and paternal homologous 
chromosomes (20). Its localization and distribution along the 
chromosomes indicates meiotic progression in PSC differenti- 
ation to germ cells (5,7,8). Ectopic overexpression of VASA in 
pluripotent human iPSCs has recently been shown to promote 
meiotic progression asjudged by positive SCP3 staining patterns 
(6). Thus, we compared expression and localization patterns of 
SCP3 protein in OSKM- and OSKMV-differentiated cells. Our 
analysis revealed that the majority of cells regardless of VASA 
mRNA addition during reprogramming stained negative for 
SCP3 protein indicating no meiotic activity (Fig. 2C and Supple- 
mentary Material, Fig. S4C). Nonetheless, clusters of differen- 
tiated cells showed a punctuated staining pattern, indicative of 
early premeiotic events (leptotene stage); moreover, cells repro- 
grammed with OSKMV displayed slightly higher percentages 
{P < 0.05) of cells positive for SCP3 staining. In addition, we 
noticed a distinct staining pattern that was most prominent in 
HI hESCs and in iPSCs reprogrammed with OSKMV consisting 
of elevated SCP3 clusters with elongated structures indicative of 
late, albeit disorganized, assembly of meiotic chromosomes in 
zygotene, pachytene or diplotene meiotic prophase I stages 
(5,7). This staining pattern was observed in ~6% of cells, a 
percentage that is significantly higher {P < 0.05) than in the 
OSKM-reprogrammed cells or the female hESC control (8 and 
10%, respectively) (Fig. 2C). Few cells in the undifferentiated 
negative female control (H9 hESC) displayed SCP3 staining 
patterns as expected. Interestingly, the second differentiated 
male control (HI hESC) line displayed the greatest number of 
cells with elongated SCP3 staining pattern comparable with 
differentiated OSKMV cells. 

Analysis of epigenetic status and global transcription profiles 

We next examined the epigenetic status and global transcription- 
al profiles of derived lines (Supplementary Material, Fig. S4D 
and E). We analyzed four imprinted genes via bisulfite 
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Figure 2. Functional and molecular studies of iPSC.HUF 1 derived with OSKM and OSKMV. (A) Gene expression analysis of markers associated with the germline 
lineage in an undifferentiated state. Three genes (PRDM 14, VASA, DPPA3) showed significantly higher expression in OSKMV-derived clones (purple bar) compared 
with OSKM (Student's Mest, mean + SD;h > 1 6 for each gene and sample *P < 0.05, **P < 0.01, < 0.001, < 0.000 l).(B) Gene expression analysis of 
clones derived with OSKM and OSKMV, iPSC.BJ line and H9 during in vitro differentiation into primordial germ cells with BMP4. Samples were isolated at day 0, 2 
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sequencing: KCNQIOTI- and PEGl/MEST-linked differen- 
tially methylated regions (DMRs), both maternally imprinted 
(2 1 ,22) as well as the H 1 9 DMR and H 1 9 promoter loci, both pa- 
ternally imprinted (23). Our analysis revealed slight but signifi- 
cant differences {P < 0.0001) in methylation status of both 
maternally and paternally imprinted loci, in particular, the ma- 
ternally imprinted KCQIOTI and the paternally imprinted 
H19 promoter region (Supplementary Material, Fig. S4E). For 
both loci, the OSKMV line showed a significantly lower percent- 
age of methylated sequences compared with the OSKM counter- 
part and the parental HUFl fibroblast line. To further test our 
findings, we also included an additional control of three inde- 
pendent pooled male sperm samples. We observed little if any 
methylation at the KCNQIOTI and PEGl/MEST loci as 
expected for maternal-specific germline methylated loci. Also, 
as expected, the paternally imprinted HI 9 DMR region demon- 
strated near-complete methylation. In contrast, to the H 1 9 DMR, 
which acquires methylation during gametogenesis, the H 1 9 pro- 
moter is further methylated following implantation (24) and was 
characterized by few methylated sequences (18%). Equal pro- 
portions of fully methylated and fully unmethylated DNA was 
detected for PEGl/MEST in both, the OSKM and OSKMV 
line {P > 0.05, test). We note that methylation patterns of 
imprinted genes are diagnostic of germ cells and that cells 
were analyzed in their undifferentiated state. 

To examine global gene expression, we performed RNAseq 
followed by differential gene expression analysis on the 
OSKM- and OSKMV-derived HUFl line together with controls 
(parental HUFl and H9). As expected, gene expression profiles 
between HUFl fibroblasts and the three pluripotent populations 
(OSKM, OSKMV, H9 hESCs) were significantly different 
(P < 0.05) (Supplementary Material, Fig. S5). Interestingly, 
HUF 1 fibroblasts and H9 hESCs expressed more genes at signifi- 
cantly different levels {P < 0.05) than HUFl fibroblasts com- 
pared with their iPSC derivatives (Supplementary Material, 
Fig. S5C). Pairwise volcano and scatter plots emphasize the 
gene expression differences between each sample (Supplemen- 
tary Material, Fig. S5A and B). We further explored the relation- 
ship between all four lines and visualized the Jensen -Shannon 
(JS) distances in a heat map and a significant genes {P < 0.05) 
overview matrix across all genes (Supplementary Material, 
Fig. S5C andD). Principal component analysis (PC A) and multi- 
dimensional scaling (MDS) revealed a very similar global gene 
expression in the OSKM- and OSKMV-derived lines that were 
significantly different (P < 0.05) from the parental fibroblasts 
(Supplementary Material, Fig. S5E and F). Derived iPSC lines 
also were distinct from the H9 hESC cell control but revealed 
a similar global gene expression phenotype (all three pluripotent 
lines with a positive second PC in PCA analysis). Hierarchical 
clustering analysis (Supplementary Material, Fig. S5G) and a 
heat map of differentially expressed genes (Supplementary Ma- 
terial, Fig. S5H) highlighted the relationship between all PSC 
lines. When we compared both derived iPSC lines with hESCs, 
we observed only a few genes, isoforms and transcription start 
sites (TSSs) with statistically significant differential expression, 
consistent with previous reports (25,26). 



Xenotransplantation directs germ cell formation in OSKM- 
and OSKMV-reprogrammed cells 

There are no previous reports of transplantation of human PSCs 
to the seminiferous environment in vivo. The use of testicular 
transplantation is common for teratoma analysis; however, as 
commonly used, teratoma formation assays rely on transplant- 
ation under the testis capsule. In contrast, here we made use of 
a xenotransplant assay in which cells from reprogrammed 
lines were directly injected into the seminiferous tubules of 
busulfan-treated immunodeficient nude mice as previously 
described (27,28). We used this assay to investigate OSKMV 
and OSKM donor cell differentiation into germ cells in vivo 
with injection directly into mouse seminiferous tubules. Specif- 
ically, we hypothesized that the inclusion of VASA during repro- 
gramming would confer a unique gerni cell competent state. As 
previously reported, transplanted human spermatogonial stem 
cells migrate to the seminiferous tubule basement membrane 
and proliferate to form chains and patches of spermatogonia 
that persist long term (28). Note that complete spennatogenesis 
is a function of evolutionary distance such that rat transplant- 
ation into mouse tubules or human transplantation into non- 
human primates is expected to yield complete spermatogenesis 
whereas human transplantation of germ cells into mouse or rat 
tubules will be limited to formation of prospermatogonia, sperm- 
atogonia and/or possibly early meiotic derivatives (28). As a 
positive control, we transplanted human fetal testicular cells 
(22-week-old tissue) into busulfan-depleted spermatogenic 
tubules of immunodeficient mice and observed clusters and 
chains of spermatogonia 2-month post-transplantation (Fig. 3A, 
left panel). Transplantation of HI hESCs (XY karyotype) and 
H9 hESCs (XX karyotype) served as controls; we hypothesized 
that H9 hESCs would result in formation of few, if any, germ 
cells in seminiferous tubules relative to the XY hESCs. Two stand- 
ard methods to examine the potential to form germ cells in trans- 
plants were pursued: Immunohistochemistry of serial sections of 
fixed tissue and whole-mount staining of testis to assess potential 
differentiation to chains of spermatogonia, in vivo. Whole mount- 
ing is only possible when engraftment occurs in the absence of 
large masses post-transplantation. For immunohistochemical 
(IHC) analysis of serial cross sections of transplanted tissues, 
we used a panel of very well-characterized germ cell markers 
(2,5,29-32). 

We injected undifferentiated OSKM and OSKMV iPSCs into 
seminiferous tubules of eight testes each to evaluate their potential 
to form PGCs, gonocytes or spermatogonia in vivo in a transplant 
assay. Depending on which factors were used during reprogram- 
ming, we obtained strikingly different results. Testes injected with 
OSKMV cells maintained their naive tissue structure 2-month 
post injection and one of eight testes was positive for a cluster 
of human cells that persisted long term (Fig. 3A, right panel and 
Supplementary Material, Table S3). In contrast, testes trans- 
planted with OSKM cells developed large internal proliferating 
cell masses with no signs of teratoma formation in all eight of 
eight transplanted testes — ^transforming the testis into an enlarged 
tissue (Fig. 4A and B and Supplementary Material, Fig. S6A). 



and4 post-BMP4/vehicle treatment. Thi-ee key markers (NANOS3, VASA, DPP A3) were upregulated in OSKMV-derived clones (Student's t-test, mean + SD; n > 
16 for each gene and sample *P < 0.05, **P < 0.01, ***p < 0.001, ****/> < 0.0001). (C) Quantification of meiotic spreads. Approximately 200-300 cells were 
counted for each line. Staining patterns were classified as negative, punctuated or elongated. 
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To examine outcomes more extensively, we then repeated the 
transplantation and used the same method of direct comparison 
of germ cell activity across all transplantations. Our replicate 
transplantations consisted of OSKMV cell injections into the sem- 
iniferous tubules of an additional six testes with subsequent IHC 
analysis of serial cross sections regardless of tissue structure. 
We confinned previous observations that indicated that OSKMV- 
derived iPSCs do not form large masses of cells but instead leave 
the mouse testis structure intact. In direct contrast, testes trans- 
planted with undifferentiated H9 (Fig. 4Av) and HI (Fig. 4Aiii) 
hESCs, developed enlarged tissues indicative of tumors similar 
to OSKM cells (Fig. 4Aii) and in addition formed teratoma-like 
structures (in H9 cells only. Fig. 4Avi) indicative of multiple 
germ layer tissue formation (see also Supplementary Material, 
Table S3). We observed that HI cells that were not localized to 
the basement membrane did not demonstrate clear differentiation 
to either germ cells or somatic cells. Instead, based on the hist- 
ology of multiple xenografts, HI and OSKM cells outside the 
tubules resembled the histology of embryonal carcinoma (EC) 
cells and yolk sac tumors (Fig. 4Avii and viii) (33). 

Next, we focused on directly comparing germ cell activity and 
extensively screened serial cross sections for a panel of well- 
established germ cell markers in all transplanted mouse testis 
(OSKM, OSKMV, HI, H9). We compared our results with 
cross sections from a 22-week-old human fetal testis sample 
which contain abundant prospermatogonia and spermatogonia. 
Our analysis depended on the identification of human donor 
cells and the ability to discriminate them from mouse cells. 
For this purpose, we used a human-specific antibody, NUMA, 
that stains the cell nucleus and validated its specificity by stain- 
ing various tissues of human and mouse origin with latter stain- 
ing negative for NUMA (Supplementary Material, Fig. S6B and 
C). Our analysis revealed germ cells inside tubules of OSKM- and 
OSKMV -transplanted testes that are of human origin (nuclear 
NUMA) and co-localize with several germ cell-specific markers 
including VASA and DAZ (Fig. 3C and D and Supplementary 
Material, Movies SI and S2). In contrast, transplanted H9 (XX 
female) cells that were detected inside tubules did not stain posi- 
tive for any germ cell markers (Fig. 3D). Note, busulfan-depleted 
seminiferous tubules bear the risk of incomplete germ cell deple- 
tion resulting in infrequent mouse germ cells present; thus, we 
carefully discriminated positively stained germ cells that are of 
human or mouse origin (Fig. 5 A) in all of our analysis and subse- 
quent quantification efforts. 

OSKMV in contrast to OSKM reveals greater germ cell 
forming potential 

In order to quantify our immunohistochemical analysis, we 
counted tubules that stained positive for all VASA/NUMA 
cells to determine the average fraction of tubules that had resid- 
ing human germ cells present (Fig. 4C). We observed 7 and 1 8% 



of tubules positive for both OSKM- and OSKMV-transplanted 
cells, respectively (Fig. 4D and E). In contrast, we observed 
over 80% of tubules filled with NUMA/VASA double-positive 
cells in the human fetal testis cross sections (notably, these are 
non-transplanted cells and naturally existing germ cells). We 
then considered only those tubules positive for NUMA/VASA 
activity and counted individual cells across all sections; results 
indicate a more than 4-fold difference between OSKM and 
OSKMV with the latter being superior (Fig. 4F). To compare 
germ cell production between OSKM and OSKMV derivations 
and our human fetal testis control, we calculated relative germ 
cell numbers by multiplying the percent occupied tubules by 
the number of cells per tubule and compared them to our positive 
control. This indicated a 5-fold difference between OSKM com- 
pared with OSKMV and the human fetal testis, which was sig- 
nificant (Fig. 4G). Analysis of HI donor-derived cells revealed 
results that were similar to OSKM-transplanted cells. As high- 
lighted from the IHC analysis, mouse testes transplanted with 
H9 cells had lower germ cell activity in our quantitative analysis, 
which, however, was not significant (Fig. 4E-G). 

In order to further verify that germ cells were induced from 
iPSCs in vivo, we extended our analysis of previously characterized 
germ cell markers (2,5,29-32) and demonstrated co-localization of 
three sets of PGC and premeiotic germ cell markers ( VASA/D AZL, 
DPPA3/UTF1 and VASA/UTFl) with NUMA (Fig. 5C, Panels 1 
and 2) for OSKM- and OSKMV-transplanted cells in conjunction 
with the positive control (Fig. 5C). In addition, we observed, in 
rare instances, cells that stained positive for GFRct 1 cells, suggest- 
ing that OSKMV iPSCs are not only superior in germ cell formation 
to OSKM iPSCs in vitro, but also in vivo. 

To evaluate whether transplanted cells differentiated to other 
lineages (such as Sertoli cells) or remained in their undifferenti- 
ated state, we co-stained for NUMA and germ cell markers with 
GATA4, a Sertoli cell marker. We observed GATA4-positive 
cells (Fig. 5B) at the basement membrane of tubules that did 
not co-stain with NUMA, indicating their murine origin; we 
note that the GATA4 cells resided in the vicinity of NUMA/ 
DAZ double-positive cells highlighting their supportive pheno- 
type in germ cell differentiation. Tubules that were entirely filled 
with large numbers of cells both intratubular and extratubular 
appeared to be unorganized, of human origin and OCT3/4 posi- 
tive, suggesting an excess of transplanted cells that retained their 
undifferentiated state (Supplementary Material, Fig. S6D). This 
is in contrast to OCT3/4 positive cells that did not co-localize 
with NUMA but instead revealed an organized staining pattern 
along the basement membrane indicating endogenous cell of 
mouse origin. We observed that a small number of the trans- 
planted iPSCs might escape during injection particularly in HI 
and H9 hESC-containing testes. Specifically, HI and H9 donor 
cells appeared to exit the seminiferous tubules and entered the 
extratubular space within the mouse testis as illustrated by posi- 
tive NUMA staining outside mouse tubules (Supplementary 



Figure 3. Transplantation of OSKM and OSKMV cells into busulfan-treated mouse testes. (A) Whole-mount analysis on transplanted human fetal testis cells and 
OSKMV cells into mouse testes. Chain (white dashed rectangle) and cluster formation (white arrows) visible in human fetal testis control cells. Transplanted 
OSKMV cells gave rise to cluster formation indicated by white dashed rectangle. Scale bar, 50 |j,m. (B) Histology cross sections of tubules inside mouse testes 
stained with co-localizing human-specific NUMA (red) and VASA (green). White an'ows indicate transplanted cells positive for VASA germ cell marker; red 
dashed arrows indicate VASA-negative transplanted donor cells. Scale bar, 50 jjim. (C) Histology cross sections of tubules of mouse testes stained with additional 
germ cell-specific marker DAZ demonstrating NUMA co-localization. Scale bar, 50 jjim. (D) HI and H9 hESC controls. HI cells are positive for VASA/NUMA 
co-staining; H9 cells are negative for genu cell marker selection. White arrows indicate transplanted cells positive for VASA gerai cell marker; red dashed aiTows 
indicate VASA-negative donor cells. Scale bar, 50 |j,m. 
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Material, Fig. S6D, middle panel). Since the cells are undifferen- 
tiated iPSCs; however, they may subsequently proliferate. 



Epigenetic analysis of donor-derived germ cells from HI 
hESCs, OSKM and OSKMV iPSCs 

In order to evaluate if germ cells derived from donor iPSCs were 
specifically undergoing epigenetic remodeling events such as 
demethylation and conversion of 5-methylcytosine (5-MC) to 
5-hydroxy-methylation (5-hMC) (34,35), we performed immu- 
nohistochemistry for 5-MC and 5-hMC in human fetal testes and 
in recipient mouse testes xenografts as previously demonstrated 
for embryo analysis (36) (Fig. 6). We observed that, in fetal 
testes, the levels of 5-MC appeared reduced in a large number 
of germ cells compared with 5-hMC levels. This was particularly 
evident in that a subset of VASA-positive germ cells had little to 
no 5-MC signal and a robust 5-hMC signal in fetal germ cells 
(Fig. 6A). We stained for NUMA/VASA in adjacent tissue sec- 
tions and observed that OSKM, OSKMV and HI hESC donor- 
derived cells exhibited a similar pattern of 5-MC and 5-hMC 
signals to fetal germ cells (Fig. 6, arrows). 5-MC levels were 
low to absent in NUMA + VASA+ donor cells indicating that 
global DNA demethylation was initiated but only partially com- 
plete owing to the maintenance of 5-hMC levels, which is 
an intermediate step in germ cell DNA demethylation. We 
also noted that NUMA-positive donor cells of OSKM and 
OSKMV xenografts that were VASA-negative and not consid- 
ered germ cells expressed very high levels of both 5-MC and 
5-hMC (Fig. 6C and D, right column). These observations 
collectively suggest an epigenetic transition from 5-MC to 
5-hMC expression in fetal germ cells which is mirrored in 
germ cells derived in vivo from undifferentiated OSKM and 
OSKMV donor cells. 



DISCUSSION 

Patient-specific iPSCs may enable genetic analysis of naturally 
occurring deletions, insertions and mutations in human germ 
cell development. However, although there has been great 
progress in differentiation of mouse PGCs from PSCs via a com- 
bination of in vitro priming and in vivo use of the niche to com- 
plete spermatogenesis and oogenesis (2-4,37-39), several 
roadblocks must be overcome in order to achieve efficient, repro- 
ducible human germ cell differentiation amenable to develop- 
mental genetic studies. Here, we tested the hypothesis that 
somatic cells could be reprogrammed to a pluripotent state 



with an additional factor, VASA, in concert with the Yamanaka 
factors. We further hypothesized that VASA would endow re- 
programmed cells with an enhanced ability to form human 
germ cells relative to cells reprogrammed with Yamanaka 
factors alone. Findings presented here provide the most viable 
strategy to date to direct human germ cell differentiation from 
PSCs for several reasons: (1) iPSCs were produced via mRNA 
reprogramming, thus insuring they are integration and xeno-free 
and do not require extensive culture (which introduces genetic 
and epigenetic mutations) to remove DNA-based reprogram- 
ming factors, (2) undifferentiated iPSCs are modified via inclu- 
sion of VASA in the reprogramming cocktail in order to enhance 
germ cell development and minimize tumorigenesis post- 
transplantation and (3) our strategy makes use of the seminifer- 
ous tubule niche in vivo in order to efficiently direct human germ 
cell formation from iPSCs, thereby potentially minimizing 
errors in erasure and establishment in sex-specific imprinted 
genes as have previously been observed in all human studies to 
date (5-8,40). 

It is notable that efforts to simply translate results, in whole, 
from mouse to humans in diverse systems from hematopoietic, 
neural, cardiac, liver and others have all demonstrated a need 
to develop a specific protocol for differentiation of human 
tissues relative to the mouse. We have found this to be the case 
to date with differentiation of germ cells from mouse and 
human PSCs. Although some methods and molecules are trans- 
ferrable, others are not and thus, we focused here on alterations to 
accommodate human germ cell differentiation. Our data reveal 
that reprogramming with inclusion of VASA can endow a 
ground state of pluripotency that is unique from OSKM and 
hESC lines. This is reflected at the undifferentiated state in 
OSKMV cells by notably lowered expression of OCT3/4 
(P0U5F1), TERT and SALL4. We suggest that the lowered 
level of these core pluripotency regulators is sufficient to main- 
tain several phenotypic features in OSKMV cells consistent with 
conventially derived OSKM iPSCs. Interestingly, ectopic 
VASA expression does not affect the global transcriptional pro- 
grams of OSKMV- and OSKM-derived iPSCs. It appears that 
ectopic VASA expression during reprogramming to iPSCs pro- 
vides a signal that downregulates genes linked to pluripotency 
but does not induce in vitro differentiation under proliferative 
cell culture conditions. Our data also reveal that the OSKMV- 
reprogramming process facilitates germ cell formation from 
these cells. Evidence for this is provided in part by significant 
upregulation of premeiotic germ cell markers NAN0S3, 
VASA and DPP A3 upon BMP4 differentiation and a greater po- 
tential to enter meiosis compared with the OSKM counterparts. 



Figure 4. Transplantation of OSKM and OSKMV followed by quantification analysis. (A) Hematoxylin and eosin staining of histology cross sections of tubules inside 
mouse testes, (i-v) Black an'ows indicate individual tubules. Yellow arrows indicate cartilage formation in H9 transplanted cells as example for teratoma formation. 
Scale bar, 100 |j,m. (vi) H9 transplanted cells formed teratomas indicated by three germ layers, (vii and viii) HI and OSKM cells formed fibrotic structures that 
resembled embryonal carcinoma (EC) cells and yolk sac tumors. (B) Weight of testis 2 months after transplantation of cells (OSKM and OSKMV cells with three 
controls: human testis, HI hESC and H9 hESC line). (C) Representative low-magnification image of histology cross section of multiple tubules stained with 
NUMA and VASA indicating that only a fraction of tubules formed VASA-positive cells. Red dashed rectangle indicates tubules with NUMA only cells along 
the basement membrane. White rectangle indicates NUMA/VASA-positive tubules. Scale bar, 50 |j,m. (D) Representative image of histology cross section of one 
single tubule. Single NUMA/VASA-positive cells were counted for quantification. Scale bar, 50 jxm. (E-G) Quantification of immunohistochemistry results of 
cross sections from OSKM and OSKMV cells compared with both controls. All serial sections were subject to counting (see Materials and Methods) Scale bar, 
80 (jLm. (E) Percentage of tubules with positive VASA/NUMA co-staining calculated (against total number of counted tubules). iPSC.OSfQVIV is also significantly 
different from H9 control (P = 0.0432) (F) For each positive tubule, VASA-positive cells that co-stained with NUMA were counted and calculated against positively 
stained tubules. iPSC. OSKMV is also significantly different from H9 control (P — 0.0397) (G) Relative germ cell forming potential calculated by multiplying fraction 
of positively stained tubules with number of VASA/NUMA co-stained cells for each sample. iPSC. OSKMV is also significantly different from H9 control (P — 
0.0417). 
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Figure 5. Genu cell foiTnation of OSKM- and OSKMV-transplanted cells in vivo. (A) Obsei'ved endogenous genu cell activity inside the basement membrane of 
transplanted mouse testes. Cells stained positive for VASA, PLZF and DAZ but do not co-localize with NUMA indicating residual germ cells after busulfan treatment. 
Scale bar, 50 jxm. (B) Detection of Sertoli cells indicated by GATA4 staining that is of mouse origin and near OSKMV-transplanted cells at the basement membrane of 
mouse testes. Scale bar, 50 [Lva. (C) Additional genn cell markers stain positive with NUMA in OSKM- and OSKMV-transplanted cells in mouse testes and co-localize 
with each other. GFRa 1 was only detected in OSKMV-transplanted cells and human fetal testis control sections. Yellow dashed rectangles indicate magnified snap- 
shots in each panel and for each sample. Each panel shows single individual protein and DAPI (left) and a zoomed out merged image (right). Scale bar, 80 (jum. 
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Figure 6. DNA demethylation in xenotransplanted hESCs and iPSCs. Human fetal testis section (22 weeks) and testis xenografts of H 1 hESCs, OSKM and OSKMV 
human iPSCs were stained for VASA, 5-metliylcytosine (5-MC) and 5-hydroxymethlcytosine (5-liMC) in NUMA+ regions. (A) Cross section of a human fetal testis 
with positive immunostaining for VASA, 5-MC and 5-hMC. Areas in dotted yellow rectangles are shown in higher magnification on the left of each panel. White 
arrows indicate germ cells with low to none 5-MC signal relative to 5-hMC. (B-D) Cross sections of mouse testes xenografts after transplantation of undifferentiated 
HI hESCs, OSKM and OSKMV cells, respectively, immunostained for NUMA, VASA, 5-MC and 5-hMC. Areas in dotted yellow rectangles are shown in higher 
magnification on the left of each panel. White arrows indicate NUMA + VASA-positive donor cells with low to none 5-MC signal relative to 5-hMC. Non-germ 
cell containing regions of the xenografts are shown in far right column for OSKM and OSKMV cells. Scale bar, 50 (ji,m. 



It is worth nothing that PGC cuhure in vitro under prohferative 
conditions can resuh in derivation of pluripotent embryonic 
germ cell lines capable of differentiating to somatic and germline 
progeny cells; moreover, several studies have provided evidence 
of extensive similarity between ESCs and early germ cells 
(11,41). In contrast, the ontogeny of mouse iPSCs versus 
mouse PGCs is a more distant one than in humans (3). Mouse 
ESCs and iPSCs are considered to have a more naive ground 
state than hESCs and therefore, an epiblast-like state is 
induced prior to PGC formation (2-4). 

Given, the close relationship between human iPSCs and PGCs 
and the more distant relationship to mouse ESCs, we tested 



whether differentiation of RNA-reprogrammed iPSCs might 
simply be directed to the germline via inj ection directly into sem- 
iniferous tubules. Our observations reveal differences in 
OSKMV versus OSKM cells in favor of a germ-cell fate. This in- 
terpretation is based on several observations: First, NUMA-I- 
OSKM and OSKMV cells in testicular cross sections integrated 
into or proximal to the basal membrane of the mouse tubules 
where spermatogonia and gonocytes normally reside in early 
fetal stages. Second, OSKM and OSKMV donor cells expressed 
key premeiotic germline markers with specific expression of the 
early spermatogonial marker, GFRal limited to OSKMV cells 
only. Third, the efficiency of derivation of germ cells was 
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Figure 7. Model of transplantation of undifferentiated iPSCs into mouse testes. A 
schematic summarizing the major findings of this study. iPSCs reprogrammed 
with OSKM and OSfCMV displayed differences in their phenotype after deriv- 
ation. When transplanted, iPSC. OSKMV cells specifically differentiate to 
PGCs inside the spermatogonial tubule niche where they contact Sertoli cells 
and were positive for a set of different stage-specific germ cell markers including 
GFRal. iPSC.OSKM cells also localized inside the mouse testis niche and 
expressed key germ cell markers but to a lower extend and frequency. In addition, 
iPSC.OSKM cells developed enlarged cell masses inside mouse testes that 
resembled embryonal carcinoma-like tissue. This is in contrast to iPSC. OSKMV 
where all transplanted mouse testes kept their naive structure. The transplantation 
strategies proposed here offer a potential avenue for fertility restoration for infer- 
tile men. 



greater with OSKMV donor cells relative to OSKM cells. 
Fourth, we detected partial and in some cases complete DNA 
demethylation of HI, OSKM and OSKMV donor-derived cells 
in a germ cell-specific manner, consistent with their progression 
through the germ cell lineage. The DNA demethylation pattern 
in donor-derived germ cells closely mirrored that observed in 
second trimester fetal gonocytes and prospermatogonia and 
offered us a hint that donor cells may be undergoing a germ cell- 
specific event. However, in order to extend our results in future 
studies, locus-specific imprinting and chromatin modifications 
will be determined and a detailed analysis of donor-derived 
cells and their equivalent fetal gem cells will be performed in 
order to pinpoint the precise stage of germ cell differentiation 
that donor cells are progressing through. Finally, and of particu- 
lar importance if we are to probe genetics of human germ cell de- 
velopment in the absence of complicating factors, cells 
reprogrammed with OSKMV formed no teratomas, tumors or 
undifferentiated cell masses in all injected sites (14 testes) in 
contrast to OSKM cells (similar to HI [XY] control cells). 

Our studies indicate that the PSC state is under control of a 
transcriptional circuitry that includes the Yamanaka reprogram- 
ming factors (42). Previous studies indicate that this transcrip- 
tional program is implemented in the context of an 'open' 
chromatin state, and it has been proposed that this state allows 
transcriptional programs to switch rapidly upon induction of dif- 
ferentiation or reprogramming ( 1 6,43). We suggest that chroma- 
tin remodelers that maintain the ESC/iPSC state, may re-open 
somatic chromatin during reprogramming thus allowing 
VASA (and the Yamanaka factors) to activate the transcriptional 
network for pluripotency and/or PGC specification. Further- 
more, we suggest that the activation of a germline differentiation 
pathway in OSKMV cells is mediated by a molecular 'switch'. 



potentially a chromatin remodeler, which detects changes in 
the cell culture environment or when cells are exposed to the 
seminiferous tubule niche. In the case of OSKMV cells, this 
'switch' guarantees a germline fate while, in OSKM cells, the 
switch may require the function of additional molecular mechan- 
isms to favor the same fate. In support of this notion, chicken 
ESCs primed with ectopic VASA (Cvh) maintained pluripo- 
tency and exhibited increased germline competence when differ- 
entiated (44). Collectively, we compare differential potential of 
OSKMV cells over OSKM cells to form PGCs in the environ- 
ment of the seminiferous tubule without the risk of tumor forma- 
tion (Fig. 7). We therefore highlight the inclusion of VASA 
during iPSC reprogramming as a defined method to induce 
germ cell fate in humans. We further speculate that other 
factors with conserved functions in germ cell specification in 
mice and humans such as the transcription factor PRDM14 (3,4) 
or the translational regulators, NAN0S3, DAZL and PUMI/2 
(5,6,45) may also be employed in mRNA-based reprogramming 
to confer higher germ cell potential on target cells. 

MATERIALS AND METHODS 

Xenotransplantation assay 

Human cell lines were transplanted into the testes of busulfan- 
treated, immunodeficient nude mice (NCr nu/nu; Taconic) as 
previously described for primate (27) and human (46) spermato- 
gonia. Brieiiy, immunodeficient nude mice were treated with a 
single dose of busulfan (40 mg/kg, Sigma) at 6 weeks of age to 
eliminate endogenous spermatogenesis. Xenotransplantation 
was then performed 5 weeks after busulfan treatment by inject- 
ing cell suspensions containing 10% trypan blue (Invitrogen) 
into the seminiferous tubules of recipient testes via cannulation 
of the efferent ducts. Approximately 7 |xl of cell suspension was 
injected per testis. 

Eight weeks after transplantation, recipient mouse seminifer- 
ous tubules were dispersed gently with Collagenase 1V{ 1 mg/ml) 
and DNase 7(1 mg/ml) in D-PBS and fixed in 4% paraformalde- 
hyde, as previously described (27,46). Clusters of human cells 
were observed on the basement membranes of recipient mouse 
seminiferous tubules by staining with a rabbit anti-primate testis 
cell primary antibody (27) and a goat anti-rabbit IgG secondary 
antibody conjugated with AlexaFluor488 (Invitrogen). All dehy- 
dration, rehydration and staining steps were carried out in 12-mm 
Transwell baskets (Coming Life Sciences) to prevent loss of sem- 
iniferous tubules. 



SUPPLEMENTARY MATERIAL 

Supplementary Material is available at HMG Online. 
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